Asymptotic behavior of the entropy of chains placed on stripes. 
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By using the transfer matrix approach, we investigate the asymptotic behavior of the entropy of 
flexible chains with M monomers each placed on stripes. In the limit of high density of monomers, 
we study the behavior of the entropy as a function of the density of monomers and the width of the 
stripe, inspired by recent analytical studies of this problem for the particular case of dimers (M = 2). 
We obtain the entropy in the asymptotic regime of high densities for chains with M = 2, ..,9 
monomers, as well as for the special case of polymers, where M — > oo, and find that the results 
show a regular behavior similar to the one found analytically for dimers. We also verify that in 
the low-density limit the mean-field expression for the entropy is followed by the results from our 
transfer matrix calculations. 

PACS numbers: 05.50.-|-q, 02.10.Ox, 02.70.-c 



I. INTRODUCTION 

In the thirties, the dimer model was introduced to 
mimic the adsorption of diatomic molecules on a crystal 
surface Later, this model was applied in the study 
of many other physical systems as ferroelectric and fer- 
romagnetic materials @, y, 0j S • The dimers can be 
modeled as chains with two basic units called monomers, 
occupying first neighbor sites of a lattice. A central ques- 
tion in the study of this model is to enumerate the num- 
ber of ways to place p dimers on lattice with N sites, 
such that the density of monomers is given hy p = 2p/N . 
The special case of full occupancy, where p — 1, was ex- 
actly solved for planar lattices, using a technique based 
on pfaffians 0, IM HI- However, the more general case 
p < 1, the so called monomer-dimer problem, is still an 
open question. Recently, an analytic solution was ob- 
tained for the case where there is a single vacancy at a 
certain site on the boundary of a two-dimensional lattice 

MM. 



On the other hand, in a previous work, two of us used 
the transfer matrix approach to calculate the entropy of 
flexible chains with M monomers, that we call M-mers, 
as a function of the density p of sites of the lattice oc- 
cupied by monomers (l^ . Since only the infinite energy 
related to the excluded volume interaction is considered, 
this is an athermal problem. The entropy was calculated 
through a numerically exact procedure for the model de- 
fined on stripes of finite widths n and infinite length, with 
periodical boundary conditions in the transverse direc- 
tion. The sequence of results for stripes of finite widths 
was then extrapolated to the two-dimensional limit by 
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means of finite-size scaling procedures. 

By using computational methods and the asymptotic 
theory of Pemantle and Wilson Kong [H, [11, , 
obtained, among other results, exact asymptotic expan- 
sions for the entropy of the dimers placed on stripes in the 
high density region. It may be appropriate to mention 
that the free energy f{p) defined by Kong is actually the 
adimensional entropy per lattice site s{p) = S{p)/ (Nks), 
which is related to the Hclmholtz free energy of the sys- 
tem through 
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where ks is the Boltzmann constant and we have set 
the constant internal energy of the system equal to zero. 
He observed that the amplitude of the first term in this 
expansion depends on the parity of the lattice width, n. 



Snip) ~ - <^(1 - p) ln(l - p), 



(2) 



where 



if : 



1 if n is odd 
i if n is even. 



The asymptotic behavior of the free energy in the low 
density limit p ^ was also studied by Kong in the 
monomer-dimer problem 

It is of interest to consider how the asymptotic results 
found by Kong are generalized if chains with more than 
two monomers are considered. In this work, using the 
transfer matrix approach to calculate the entropy for fi- 
nite stripes with width n, we obtain estimates for the 
amplitude ip not only for the dimer case, but also for 
larger chains M — 3,4, ...,9. We find that in general 
expression ([2]) describes well the behavior of s(p) in the 
high density limit, but the number of different values of 
the amplitude, as well as their numerical values, change 
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as chains with different molecular weights M are consid- 
ered. Also, we study the low density limit for the same 
values of M mentioned before. 

This paper is organized as follows. In section |TT] we 
present the expressions used to estimate the amplitudes 
Lp and we discuss how the transfer matrix for the problem 
is defined and obtained. The numerical results may be 
found in section UIIl and the conclusions are presented in 
section IIVI 



II. USEFUL EXPRESSIONS, DEFINITION OF 
THE TRANSFER MATRIX AND ITS 
CONSTRUCTION 



thus, using the equation ^ we obtain 

j \nz,,{p')dp' = An{l^ p)-^n{l- p)H^- p). (9) 

Deriving this last equation with respect to p, we have 

lnz„ = Cn - ^Pn ln(l - p), (10) 

where C„ = An — (fin- This expression was useful to 
obtain evidences that for all cases we studied the asymp- 
totic behavior supposed in equation ([5]) is valid, allowing 
us to estimate the amplitudes ipn 



Although it is rather natural to study the system in the 
canonical ensemble, where the number of monomers on 
the lattice is fixed, to apply the transfer matrix technique 
it is convenient to allow this number to fluctuate. We 
thus define the grand-canonical partition function as 



5(z)=^z''*^r(M,7V,p), 
p 



(3) 



where z is the activity of a monomer and T(M, N,p) is the 
number of ways to place p chains with M monomers each 
on the lattice with N sites. The density of monomers may 
be written as 



p{z) = z-^{z), 



(4) 



where ipi^) is the thermodynamic potential per lattice 
site, defined as 



(z) = lim ^ln^(z). 



(5) 



In the thermodynamic limit we may use a Legendre 
transformation to rewrite the potential as 



(j){z) ^ maxjplnz -|- s{p)}, 
p 



(6) 



which implies that 



ds 

— = - Inz, 
dp 



and therefore the entropy will be given by 



s{p) = )dp , 

Jo 



(7) 



A. Transfer matrix 

To build the transfer matrix for this problem we define 
a stripe of width n on the square lattice in the plane 
(x, y), so that 1 < x < n and — oo < y < oo. The position 
of a site may be defined by {x,y), where I < x < n and 
y are integer numbers. Periodic boundary conditions are 
assumed in both directions. Generalizing the prescription 
due to Derrida |17| for the transfer matrix of an infinite 
chain placed in stripes, we define the state of a set of n 
vertical bonds of the lattice connecting the sites at ?/o — 1 
to the sites at yo specifying the number of monomers 
already connected to this bond located on sites with y < 
yo (in the range [0, M — 1]) and the pairs of bonds which 
are connected to each other by a chain whose monomers 
are all located at sites with y < yo. This last information 
is essential to prevent the presence of rings in the allowed 
configurations. With the information above about the 
configuration of the vertical bonds located between yo ~ 1 
and yo, we may find all possible configurations of the 
vertical bonds between yo and yo + li thus defining a 
transfer matrix. Actually, it is not difficult to develop an 
algorithm for the steps involved in this procedure, which 
allows us to obtain the elements of the transfer matrix 
exactly. Restriction in memory and computer time set 
an upper limit to the values of n and M we are able to 
handle, since the size of the transfer matrix grows very 
fast as they increase. More details about this procedure 
may be found in the previous paper [l^ . 

Once the transfer matrix T is obtained, we may find 
the entropy in the thermodynamic limit using the largest 
eigenvalue of this matrix A. The grand-canonical par- 
tition function is related to the transfer matrix by the 
expression 



S(z) = rr(TO, 



(11) 



with s(0) = 0. 

Now, we may attempt to generalize the high density 
expansion for dimers and suppose that the behavior of 
the entropy in this region is given by the expression 



Snip) - ■Sn(l) 
(1-p) 



An - (fin ln(l - p). 



(8) 



where £ is the length of the stripe and we adopt periodical 
boundary contitions in the longitudinal direction as well. 
The density of monomers in the thermodynamic limit 
£ oo, p{z), will then be given by 

p{z)= lim |-^lnS(z) = --^InA. (12) 

' N^oo N dz ^ ' ndz ^ ' 
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Thus, using the expressions (|TC)|) and (|12p and supposing 
that the behavior in the high density hmit of the entropy 
for chains with M monomers is given by the relation 
([5]), we may estimate the amphtudes f for the set the 
molecular weights M widths of the stripe n we were able 
to handle. 



III. NUMERICAL RESULTS 

A roughly exponential growth of the size of the trans- 
fer matrix with both the molecular weight M and the 
width of the stripe n prevents us from obtaining results 
for larger chains or stripes. In this paper we show results 
for chains with molecular weight ranging between 2 and 
9 and widths of the stripes ranging between 2 and 12 
(dimers case) and 2 (M = 9 case). First, we analyze the 
high density limit behavior of the entropy for the molec- 
ular weights and widths of the stripes mentioned. Later, 
we will turn our attention to the low density limit for 
the same cases. For each pair of values for n and A/, 
we obtain the elements of the transfer matrix using the 
algorithm mentioned above in a numerically exact fash- 
ion. Then, we take advantage of the C„ symmetry of 
the states and use the power method to find the largest 
eigenvalue of the matrix A and its derivative with respect 
to z. 



A. High density limit 

We start with our results for dimers, which are in agree- 
ment with the values obtained by Kong [l^, [l^, with 
(ys = 1, if n is odd and ip — 1/2, if n is even, as is shown 
in figure [1] 
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FIG. 1: Results for the amplitude ifi for dimers (M = 2). 
At left panel the widths of the stripes are even values ranging 
between n — 2— 12, while at right panel results for odd values 
of the widths between n = 3 — 11 are displayed. The dashed 
lines indicates slopes equal to the known value of ip in each 



For chains with M = 3 (trimers) , our results also lead 
to two values of the amplitude: ip = 1/3 for widths that 



are multiples of 3 and (p = 1 otherwise. Our results for 
tetramers M — 4 show three values for the amplitude (p, 

{1 if n is odd 
i if n is a multiple of 2, but not of 4 
J otherwise. 

Finally, in the case M = 5 (pentamers), we found (p = 
1/5, if n is a multiple of 5 and ip — 1, in all other cases. 
Unfortunately, in this case we only were able to consider 
widths up to n = 5, due to computational limitations 
caused by the fast growth of the transfer matrix with n. 
All these results are summarized in figure [2l 



In our calculations, numerical errors set an upper limit 
to the densities we may consider. As we approach the 
full occupancy limit p — > 1, the values of the activity z 
become very large and this makes the numerical errors 
in the calculations grow. Finally, we extend our analysis 
to the polymer case (M — * oo). As was shown in [12 |. 
the entropy for this case shows different finite size scaling 
corrections for even and odd widths. Nevertheless, the 
data show a single amplitude </? = 1 in this limit, as may 
be seen in figure [H 



We may conclude that the results shown above sug- 
gest that the dependence of the entropy with the width 
of the stripe obtained by Kong for the dimer model may 
be extended for larger chains, including more possibili- 
ties for the amplitude (p. Apparently, the values of the 
amplitude are related to the split of the leading eigen- 
values of the transfer matrix into subsets with different 
finite size scaling behaviors found in [12]. These splits 
seem to be related to frustration effects in the limit of 
full occupancy. An empirical rule that we found to label 
these subsets is the determination of an integer a which 
minimizes the relation 



(13) 



where k is the smallest integer for which a solution is 
found. The widths associated to the same value of a 
share the same estimate of the amplitude, as may be 
seen in table HI Another way to look at relation ([TH)) is to 
rewrite it as a/k = M/n, so that a is the numerator of 
the fraction M/n after it is simplified. We notice that a 
may be interpreted as the length of the smallest rectangle 
of width n which may be totally filled by chains with M 
monomers each. The number of chains we can place in 
this rectangle is k. We notice that if M is prime, we 
find a — M for all widths n which are not multiples of 
M, while a = 1 if n is a multiple of M. In general, the 
number of different values for a is equal to the number 
of divisors of M, including 1 and M itself. 

Our results suggest that the amplitude (p is related to 
the integer a by the simple relation, 



a 



(14) 
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FIG. 3: Polymer case (M oo), in which only one amplitude 
if — 1 was found for all widths studied (n = 2, 3, . . . , 9). 



TABLE I; Values for the integer a, which satisfies the relation 
l|13p for some widths and molecular weights. 
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FIG. 2: Plots which lead to the estimated values of the am- 
plitude ip for M = 3, 4 and M = 5. Top: Results for trimers 
(M = 3), with widths ranging between n = 2 and n — 10. In 
this case for widths that are multiples of 3 we find Lp — 1/3 
for the all other widths the results lead to = 1. Only the 
extreme cases in each case are labeled. Middle: Case M = 4, 
for widths in the range n = 2 — 8. In this case three values of 
the amplitude are found. For widths which are even but not 
multiples of 4, ip — 1/2. If n is odd, ip = 1 and <^ = 1/4 in all 
other cases. Bottom: Data for A4 = 5, with n = 2 — 5. Again, 
values for ip are found. If n is a multiple of 5, then p — 1/5 
and 1/5=1 for the other cases. The dashed lines are drawn 
with the conjectured slope p in each case. 



which is obeyed by all numerical results presented here. 
However, we have no physical or mathematical argument 
to justify this hypothesis. Nevertheless, the numerical 
evidences support this relation. If we suppose that the 



polymer limit M ^ cx3 is approached by a sequence of 
prime values of M, the we would conclude that a = M 
in this limit, which is consistent with the observed am- 
plitude if — I. 

In order to test the conjecture above for the amplitude 
of the high density limit asymptotic form of the entropy, 
we extended the transfer matrix calculations to higher 
values of the molecular weight M, although with increas- 
ing values of M we are restricted to decreasing maximum 
widths n. As may be seen in figureUl all results for chains 
with molecular weights M = 6, 7, 8, 9 are consistent with 
the conjecture for the amplitudes (p. There are cases, like 
M = 6, n = 3, where numerical errors prevented us from 
obtaining results at densities high enough to observe the 
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asymptotic behavior. An example where the asymptotic 
regime also was reached only at very high densities may 
be seen in results for M = 3, n = 9 in figured As a con- 
sequence of these limitations, only some of the possible 
values for the amplitude Lp are observed in the results for 
chains with larger molecular weight presented in figure ID 
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FIG. 4: Results for the entropy in the high density limit for 
chains with molecular weight M between 6 and 9, in increas- 
ing values from left to right. The observed values of the am- 
plitudes are consistent with the conjecture made in the text. 



asymptotic behavior 




M 
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In order to test expression (|17p in the low density limit, 
we built curves for different molecular weights and widths 
of stripes. Our results, as shown in the figure [5l are 
consistent with the conjecture for all the chains analyzed 
here, since curves of Ms{p)/p as a function of \n{2p/M) 
are linear with a slope equal to 1 in all cases. 
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FIG. 5: Results for the low density limit with A/ varying 
between M = 2 and M — 6. In all cases, independently 
of the width analyzed, we observe linear relations with slope 
1/ = 1 (indicated by the dashed line), satisfying the mean-field 
prediction. 



B. Low density limit 



In the low density limit p ^ 0, Kong [T^, [l^ obtained 
from an asymptotic expansion the following behavior for 
the entropy of dimers placed in a stripe with width n, 

s„(p)^-^lnp + 0(p). (15) 

The expression above resembles the one predicted by the 
mean- field approximation [isl . [l9l | , which is 

sip) = _(l-p)ln(l-p)-£:ln(^H^) + 

+ (^1-^) (lng-1), (16) 

where q is the coordination number. For small values 
of the density p, the leading contribution comes from 
the second term, which for dimers is identical to the one 
obtained by Kong. 

It is then rather natural to conjecture that for other 
values of M the leading contribution to the entropy is 
the one predicted by the mean field approximation, since 
for small densities the interchain interactions may be ne- 
glected. Thus, we conjecture that for all values of the 
molecular weights in the low density limit we have the 



IV. CONCLUSIONS 

In this work, we study the asymptotic behavior of 
the entropy for chains placed on stripes using the trans- 
fer matrix approach. Generalizing the results by Kong 
[flSjIIB] for the case of dimers in the high density limit, 
we conjecture similar asymptotic forms for the entropy 
of chains with larger molecular weights. We propose an 
empirical rule to determine the number of different am- 
plitudes of the asymptotic behavior close to the full occu- 
pancy limit, as well as to find the values of these ampli- 
tudes. Unfortunately, computational limitations prevent 
us from studying larger chains and larger widths. How- 
ever, it seems reasonable expect that our results apply 
also to other cases. An analytic approach similar to the 
one made by Kong fisl . U^, using the Pemantle- Wilson 
asymptotic theory [l^ would be very useful to support 
our conclusions. However, it may not be easy to carry 
out this task for chains with molecular weight larger than 
M = 2. 

Another interesting result is that in the low density 
limit the free energy is well fitted by the mean-field ap- 
proximation p^ . [l9j , for any width n for all the values 
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of M wc were able to consider. Again, it would be inter- 
esting to obtain an exact result in this limit in order to 
verify this conclusion. 
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